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Abstract 

Motivated by the successful synthesis of several molecular quantum spin rings we are investigating whether such 
systems can host magnetic solitary waves. The small size of these spin systems forbids the application of a classical or 
continuum limit. We therefore investigate whether the time-dependent Schrodinger equation itself permits solitary 
waves. Example solutions are obtained via complete diagonalization of the underlying Heisenberg Hamiltonian. 
© 2008 Elsevier B.V. All rights reserved. 
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1. Introduction 

Magnetic solitons are detected in many magnetic 
systems due to their special influence on magnetic ob- 
servables [1,2,3,4,5]. From a theoretical point of view 
magnetic solitons are solutions of non-linear differen- 
tial equations, e. g. of the cubic Schrodinger equa- 
tion [6,7,8], which result from classical approximations 
of the respective quantum spin problem. The cubic 
Schrodinger equation for instance is obtained if the 
spins are replaced by a classical spin density [6] . 

Due to the success of coordination chemistry one 
can nowadays realize finite size quantum spin rings in 
the form of magnetic ring molecules. Such wheels - 
Fee, Feio, and Crg rings as the most prominent exam- 
ples - are almost perfect Heisenberg spin rings with a 
single isotropic antiferromagnetic exchange parameter 
and weak uniaxial anisotropy [9,10,11,12]. 

The aim of the present article is to discuss whether 
solitary waves can exist on such spin rings and if they 
do, how they look like. The finite size and the resulting 
discreteness of the energy spectrum forbid any classical 
or continuum limit. We therefore investigate, whether 
the ordinary linear time-dependent Schrodinger equa- 
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tion allows for solitary waves. Solitary excitations in 
quantum spin chains of Ising or sine-Gordon type have 
already been discussed [7]. Nevertheless, such soliton 
solutions are approximate in the sense that a kind of 
Holstein-Primakov series expansion is applied, and the 
results are accurate only for anisotropies large com- 
pared to the spin-spin coupling [7]. This article deals 
with antiferromagnetic Heisenberg chains without 
anisotropy, where such derivations cannot be applied. 
Starting from the time-dependent Schrodinger equa- 
tion therefore automatically addresses the questions 
how soliton solutions might be approached starting 
from a full quantum treatment - a question that to 
the best of our knowledge is not yet answered [13] . 

In order to apply the concept of solitons or solitary 
waves to the linear Schrodinger equation some redef- 
initions are necessary. The first redefinition concerns 
the term soliton itself. It is mostly used for domain- 
wall like solitons which are of topological character. It 
is also applied to localized deviations of the magneti- 
zation or energy distribution (envelope solitons), here 
one distinguishes between bright and dark solitons. We 
will generalize this class of objects and speak only of 
solitary waves in the following. We call a state ] '^s ) 
solitary wave if there exists a time r for which the time 
evolution equals (up to a global phase) the shift by one 
site on the spin ring. This means that solitary waves 
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travel with permanent shape. The property that two 
solitons scatter into soliton states cannot be used as a 

definition in the context of the Schrodinger equation 
since it is trivially fulfilled for a linear differential equa- 
tion. 

A nontrivial question concerns useful observables in 
order to visualize solitary waves | ^'s ). The expecta- 
tion value of the local operator Sz(i), which reflects a 
local magnetization is used classically and meaning- 
ful also as a quantum mechanical expectation value. 
Then nontrivial, i. e. non-constant magnetization dis- 
tributions { I Sz{i) \ *s ) arise only for those solitary 
waves which possess components with a common to- 
tal magnetic quantum number M, because matrix 
elements of s^z{i) between states of different total 
magnetic quantum number vanish. The energy density 
which in classical calculations is also used to picture 
solitons, could quantum mechanically be defined as 
s(i) • s(i + 1) starting from a Hciscnbcrg Hamiltonian, 
compare Eq. (1). It turns out that this observable is 
not so useful since it is featureless in most cases be- 
cause the off-diagonal elements for energy eigenstates 
belonging to different total spin S are zero. 

The article is organized as follows. In Sec. 2 we 
shortly introduce the Heisenberg model and the used 
notation, while in Sec. 3 solitary waves are defined. 
We will then discuss the construction and stability of 
(approximate) solitary waves in Sec. 4. A number of 
example spin systems will be investigated in Sec. 5. 
The article closes with conclusions in Sec. 6. 

2. Heisenberg model 

The Hamilton operator of the Heisenberg model with 
autifcrromagnetic, isotropic nearest neighbor interac- 
tion between spins of equal spin quantum number s is 
given by 

JV 

H^^Y, s{i)-s{i + l), (1) 

i=l 

H is invariant under cyclic shifts generated by the shift 
operator T. T is defined by its action on the product 
basis I m ) 

T I mi, ... , TUN ) = I mjv, mi, ... , mjv-i ) , (2) 

where the product basis is constructed from single- 
particle eigenstates of all Sz(i) 

Sz{i) I mi, . . . ,mN ) = m, I mi, . . . ,mjv ) • (3) 

The shift quantum number k = 0, . . . , N — 1 modulo 
N labels the eigenvalues of T which are the A^-th roots 
of unity 



2 

^ = exp|-i^| . (4) 

k is related to the "crystal momentum" viap = 2'nk/N . 

Altogether H, T, the square , and the z- 
component Sz of the total spin are four commuting 
operators. 

3. Solitary waves 

Wo call a state | s ) solitary wave if there exists 
a time t for which the time evolution equals (up to 
a global phase) the shift by one site on the spin ring 
either to the left or to the right, i. e. 

l/(r) |*.)=e-**°r±M*.)- (5) 

Decomposing | ^ s ) into simultaneous eigenstates 
I ) of fl" and T, 

1*,)=^ C.I*.), (6) 
yields the following relation 

^ e-*^c.|*.)=e-'*« ^ e^*^cj*^)(.7) 

The set Is contains the indices of those eigenstates 
which contribute to the solitary wave. Therefore, 

:^=±^+27rm^ + $o V/. € 7. (8) 
with rrin e Z . 

Equation (8) means nothing else than that solitary 
waves arc formed from such simultaneous eigenstates 
I 4". ) oi H and T which fulfill a (generalized) linear 
dispersion relation. 

Since our definition (8) is rather general it comprises 
several solutions. Some of them are rather trivial waves 
whereas others indeed do possess soliton character: 

- Single eigenstates | *.) of the Hamiltonian fulfill the 
definition, but they are of course stationary and they 
possess a constant magnetization distribution. One 
would not call these states solitary waves or solitons. 

- Superpositions of two eigenstates | ) and | ) 
with different shift quantum numbers fe^ and ku are 
also solutions of definition (8) since two points are 
always on a line in the S-fe-pIane. It is clear that such 
a state cannot be well localized because of the very 
limited number of momentum components. Never- 
theless, these states move around the spin ring with 
permanent shape. 

An example of such solitary waves are superpo- 
sitions consisting of ground state and first excited 
state which have already been investigated under a 



Jiirgen Schnack & Pavlo Shchelokovskyy / Journal of Magnetism and Magnetic Materials (2008) 1-0 



3 



different focus [14]. The characteristic time t to move 
by one site is related to the frequency of the coherent 
spin quantum dynamics discussed for these special 
superpositions [14,15]. 

The existence of more than two eigenstates with lin- 
ear dispersion relation cannot be proven for Heisen- 
berg models in contrast to the Ising model [7]. In 
fact it is very likely that such a combination does not 
exist at all for a small system. Therefore, superpo- 
sitions of three or more eigenstates will - hopefully 
only slightly - deviate from the linear dispersion re- 
lation and therefore describe (slowly) decaying soli- 
tary waves. This will be explained in detail below. 
The character of the solitary waves assembled ac- 
cording to definition (8) can be different. As section 5 
will show some of them are localized magnetization 
distributions (envelope solitons) which could be de- 
noted as bright or dark solitary waves others are of 
topological nature as for example in the case of su- 
perpositions of the ground state and the first excited 
state on an odd membered spin ring. Nevertheless, 
the distinction is not sharp. As shown in Ref. [16] 
envelope solitons can be of topological nature when 
looked at from a different perspective. 



4. Construction and stability of solitary waves 

When constructing solitary waves we focus on two 
principles: 

- Energetically low-lying solitons can be formed by su- 
perimposing the ground state with low-lying excited 
states. For this purpose it is advantageous that for 
antiferromagnetic rings the lowest energy eigenval- 
ues depend on momentum approximately via a sine 
function. Therefore, for large systems, where the sine 
function can be approximated by its argument, an 
almost linear dependence between and is ex- 
pected. 

- General solitons can be formed by superimposing 

states with arbitrary pairs {k^,E^) as long as they 
fulfill (8). For this purpose two states are selected 
which define a straight line in _E — fc-space. Then a 
search algorithm is used to find other states whose 
energy and momentum quantum numbers constitute 
a point on (or in close vicinity) of the line. The search 
extends not only to the first Brilluoin zone, but also 
to higher zones. 

It is clear that in a finite system the linear disper- 
sion relation (8) might only be fulfilled approximately 
which leads to decaying solitary waves. That the time 
evolution must be recurrent in a Hilbert space of finite 
dimension is only of principle interest since the recur- 
rence time would be very long. 

A good measure of the stability of a solitary wave is 
the overlap of the time-evolved state and the shifted 



one. For perfect stability these two states coincide. 
Let's assume that the wave travels in the same direc- 
tion as the shift operator acts, then we define the star 
bility measure as 



r]ir) = {^s\T-'U{r) j > 



(9) 



77(r) is a complex number, and for dispersionless move- 
ment it's absolute value equals one. Decomposing the 
approximate solitary wave | ) into simultaneous 
eigenstates \ of H and T yields 



7(r)=E 



|2 - 

\c„\ e 



(10) 



where S^, is the energy mismatch of energy level v. In 
order to consider a motion across n sites the exponent 
has to be multiplied by n. 

As a simple example we want to consider an approx- 
imate solitary wave which consists of three eigenstates 
as for instance in the example of Fig. 4. The first two 
levels define the linear dependence, and the third level 
deviates by 63 from that dependence. Then, taking the 
normalization condition into account, the overlap turns 
out to be 



r,(r) = l-|c3|' + |c3|V*^ 



(11) 



Thus we face a decaying overlap which later oscillates 
back to one. For the specific example of Fig. 4 we can 
estimate how much time this will take. The character- 
istic time T is given by 



h _ AE N 
T ^ Afc 27r 



(12) 



which is approximately 3 in the given units, compare 
Fig. 1. The energy mismatch is of the order of 10~^. 

In order to have a deviation of the exponential in (11) 
of 1 % from unity the argument has to grow up to 0.1, 
which means that the solitary wave has to move across 
about 30 sites. Or in other words, the approximate 
three state solitary wave of Fig. 4 has to circumvent the 
spin ring three times in order to produce a deviation 
of the overlap with the original state of about one per- 
cent. This means that such approximate solitary waves 
are rather stable. Nevertheless, with increasing energy 
mismatch the situation will be worse. 



5. Special solutions for spin rings 

In this section we discuss some typical examples for 
magnetic solitary waves on finite spin rings. We have 
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chosen examples which are as large as numerically pos- 
sible, i.e. N = 11, 12, 14 and s = 1/2, because a too 
small size yields only a very limited number of eigen- 
states with diflterent shift (i.e. momentum) quantum 
numbers that can be superimposed to form the solitary 
wave. Then, the solitary waves obviously cannot be of 
such perfect shape and localization as known from con- 
tinuous models. Nevertheless, our general ideas are ap- 
plicable as well to systems with spin quantum numbers 
s > 1/2, which in future investigations might be ac- 
cessible by means of Density Matrix Renormalization 
Group (DMRG) techniques [17]. 
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Fig. 1. Energy spectrum of a antiferromagnetic spin ring 
with N = 12 and s = 1/2: Tlie solid as well as the dashed 
line connect states which arc superimposed to form a soli- 
tary wave. The numbers label these states. 

Figure 1 shows a typical spectrum for an antifer- 
romagnetically coupled Heisenberg ring. The ground 
state of such a ring is a singlet with fc = for the spe- 
cial case oi N = 12. The first excited state is a triplet 
with k = Q, compare the theorems of Licb, Schultz, and 
Mattis [18,19,20] for quantum numbers of ring systems 
with an even number of sites and Ref. [21] for quantum 
numbers of ring systems with an odd number of sites. 
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Fig. 2. Solitary wave for N = 12 and s = 1/2 consisting of 
ground state and first excited state with M = (crosses). 
The curve is drawn as a guide for the eye. In the course of 

time such a wave moves around the spin ring without dis- 
persion. The local magnetization distribution is the quan- 
tum expression of the classical Neel state. 



Two basic types of low-energy solitary waves can be 
formed. The first one is a superposition of the ground 
state and the triplet state with M = (dashed line in 
Fig. 1). If added up the result is a sinusoidal magne- 
tization function - sometimes also called spin-density 
wave - which is depicted in Fig. 2. A subtraction leads 
to an inverted local magnetization. In the course of 
time these distributions just move along the spin ring 
without dispersion. This dynamics is equivalent to an 
oscillation of the local magnetization, which was in- 
tensely discussed under the aspect of tunneling of the 
Neel vector [14,22,15]. 




i 



Fig. 3. Solitary wave for TV = 11 and s = 1 consisting of 
ground state (fc = 0) and first excited state (fc = 5) with 
M = (crosses). The curve is drawn as a guide for the 
eye. The local magnetization distribution is the quantum 
expression of a topological solitary wave, where the Neel 
up-down sequence is broken and continued with a displace- 
ment of one site. This state propagates without dispersion. 

Although these states look rather unspectacular the 
same procedure applied to rings of an odd number of 
spins leads to topological solitary waves. The reason 
is that due the periodic boundary condition the anti- 
ferromagnetic order is frustrated in such cases. Figure 
3 shows the case of eleven spins with s = 1. Here the 
ground state has 5 = and fc = and the first excited 
one, which is sixfold degenerate, has 5* = 1 and k = 5 
as well as A: = 6. The superposition which is depicted 
in Fig. 3 consists of the non-degenerate ground state 
and the (M = 0, fe = 5)-component of the first excited 
state. One realizes that this state describes a displace- 
ment of the classical Noel order by one site. As sketched 
by the arrows the Neel order starts at the l.h.s. as a 
sequence where up spins are located at even sites and 
down spins at odd sites. This sequence is separated by 
a "domain wall region" from the sequence at the r.li.s. 
where up spins are located at odd sites and down spins 
at even sites. Such a behavior of the local magnetiza- 
tion can equally well be pictured as a Mobius strip, 
compare Ref. [23], although the system discussed in 
that reference does not possess translational symme- 
try. 
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Fig. 4. Several solitary waves depending on contributing 
eigenstates of H and T with M = 0. All states contribute 
with the same weight in this presentation. All states with 
more than two components disperse slowly according to 
(12). Wc estimate that the overlap with the initial state 
decreases by about 1 % after two cycles around the ring. 

The second basic type of low-energy solitary waves 
is built of states along the solid line in Fig. 1. Adding 
up two to five components with the same amplitude 
yields the magnetization distributions shown in Fig. 4. 
These solitary waves are similar to two domain walls 
on a larger ring system which separate a spin-up region 
from a spin-down region (Two domain walls are needed 
because the ring has an even number of sites.). For 
a small system such as the investigated one the size 
of the two domain walls is similar to the size of the 
whole system. The solitary waves in Fig. 4 with more 
than two components are not stable because they fulfill 
the linear dispersion relation (8) only approximately, 
compare Sec. 4. 
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Fig. 5. Solitary wave built of six components for N = 14, 
s = 1/2, and M = 2 [24]. The deviation of the energy levels 
from a strict linear behavior is less than 10~*, which means 
that this state disperses very slowly, i.e., it can perform 
many cycles around the ring without noticable deviation 
from the initial shape. 

The above discussed low-energy solitary waves have 
a magnetic quantum number M = and are acces- 
sible iu low-ternpcrature experiments. Of course, soli- 
tary waves can be constructed also for larger magnetic 
quantum numbers M although it might not be easy to 



excite or detect them. Figure 5 presents an example of 
a rather well localized solitary wave in a subspace with 
non-zero magnetic quantum number [24] . 



6. Conclusions 

Summarizing, we have presented a framework to dis- 
cuss solutions of the time-dependent Schrodinger equa- 
tion which are of permanent shape, i.e. travel without 
dispersion. Besides solutions which have a close aifinity 
to solitary waves our rather general definition (8) also 
comprises solutions which arc simpler waves, which 
nevertheless propagate without dispersion. In general 
one can say, that due to the two reasons that the in- 
vestigated systems are small and that the Heisenberg 
Hamiltonian (1) does not possess any anisotropy the 
observed solitary waves are rather broad. 

The special case of a superposition of the ground 
state and the first excited state of a spin ring should 
be well observable at low temperatures for instance by 
Nuclear Magnetic Resonance (NMR) . For rings with an 
even number of sites this is currently aimed at [14,15]. 
Rings with an odd number of sites are difficult to pro- 
duce due to steric arrangement problems of the ex- 
tended ligands. First attempts resulted in rings where 
one paramagnetic ion is different from the others [23]. 
The prospects of interesting features due to frustration 
nevertheless fuel future efforts to synthesize odd rings. 

Another important aspect is the role of anisotropic 
terms in the Hamiltonian. Such terms are relevant for 
the magnetic properties of several magnetic molecules, 
e.g. for ferric wheels comprising six iron ions of spin 
s = 5/2 [25,26,22]. Easy axis anisotropy would for in- 
stance amplify the magnetization oscillations, which 
arise when ground state and first excited state are su- 
perimposed. Figure 6 shows in its bottom panel such 
oscillations again, this time for a hexanuclear iron ring. 
One notices that the local magnetization reaches at 
most an amplitude of 1.6 due to strong quantum fiuc- 
tuations in the Heisenberg model. Increasing easy axis 
anisotropy would result in more Ising-likc behavior 
with increased local magnetizations. It will be the sub- 
ject of further studies how that in general would influ- 
ence solitary waves, but it is already clear that under 
such conditions solitary waves would be more localized 
[27]. 
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Fig. 6. Top panel: Low-lying energy spectrum of a antifer- 
romagnetic spin ring with TV = 6 and s = 5/2: The dashed 
line connects the ground state and the first excited state. 
Bottom panel: Solitary wave for N = 6 and s = 5/2 consist- 
ing of ground state (fc = 3) and first excited state {k = 0) 
with M = (crosses). The curve is drawn as a guide for 
the eye. 
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